Using time reversal symmetry for sensitive incoherent matter-wave Sagnac 

interferometry 
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We present a theory of the transmission of incoherent guided matter-waves through Sagnac inter- 
ferometers. Interferometer configurations with only one input and one output port have a property 
similar to the phase rigidity observed in the transmission through Aharonov-Bohm interferometers 
in coherent mesoscopic electronics. This property is connected to the existence of counterpropa- 
gating paths of equal length and enables the operation of such matter-wave interferometers with 
incoherent sources. High finesse interferometers of this kind have a rotation sensitivity inversely 
proportional to the square root of the finesse. 
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The Sagnac effect, occurring in a wave propagating 
through a closed rotating ring, induces a phase shift 
proportional to the angular frequency fl of this rota- 
tion and the area A of the ring 1]. For light waves 
with frequency u> this phase shift is (fright = (2LuA/c 2 )fl, 
where c is the speed of light. For de-Broglie waves of 
non-relativistic particles with mass m, thephase shift is 
^matter = (2mA/h)Q = (mc 2 1 'hoj) </>ii g ht @, 13], so that 
the rotation sensitivity of a matter-wave interferometer 
('ifm') is potentially better by a factor of ~ 10 11 . First 
experimental attempts demonstrated rotation sensitivi- 
ties comparable to or even better than those of optical 
Sagnac ifm's 0, 0, 0]- However, these "one-pass" ifm's 
are limited by their small effective area and the relatively 
low flux available from coherent matter wave sources. 
Light ifm's with "multi-pass" configurations ("high fi- 
nesse"), such as ring laser gyros, may be used to increase 
the effective area and achieve good rotation sensitivity. 
Recently waveguide ring structures for cold atoms were 
demonstrated [3, S], opening the door for the develop- 
ment of " multi-pass" guided atomic Sagnac interferome- 
try. However, the small de-Broglie wavelength and short 
coherence length (of the order of 1/im) of matter waves 
makes these ifm's very sensitive to imperfections in the 
guiding potential, leading to effective path length fluctu- 
ations. 

In this Letter we use the analogy [§] between the 
Sagnac effect for massive neutral particles and the 
Aharonov-Bohm (AB) effect in coherent electron trans- 
mission through mesoscopic rings [1(3]. AB ifm's with 
only two ports connected to the ring show the effect of 
" phase rigidity" of the transmission pattern as a function 
of the magnetic flux $, when an effective path length 
difference between the ifm arms is introduced (IH 
The locking of the transmission pattern at <I> = stems 
from the time reversal invariance and current conserva- 
tion properties of any system with two ports, which im- 
plies that transmissivity is invariant to magnetic field in- 
version <I> — * — $. We show that this property, when 
applied to matter wave Sagnac ifm's, leads to their ro- 
bustness under effective path length differences and en- 
ables their operation in a high finesse configuration and 



with incoherent sources, which are available with higher 
particle flux. Ifm's of this kind will enable increasing the 
rotation sensitivity and their miniaturization onto atom 
chips [l3| while maintaining the required sensitivity. 

In general, a Sagnac ifm consists of a loop and one 
or more junctions, each consisting of one or more beam 
splitters (BS) connecting the loop to input and output 
channels. Here we consider a one dimensional (Id) model 
in which the particles are guided in a single transverse 
mode of the waveguide. A linear junction with n ports is 
represented by a n x n unitary scattering matrix S con- 
necting the output amplitudes to the input amplitudes at 
the ports. As demonstrated in fig.l we denote the indices 
of the ports of the input junction connected to the ifm 
loop by a and (3 and the ifm input port by i. The cor- 
responding input and output amplitudes at the internal 
ports a and /3 will be denoted by a± and b± respec- 
tively and the ifm input amplitude by aj„ (normalized to 
&in = !)• The relation between the amplitudes is then 
given by 



S(3a Spfi 



(i) 



If the system is linear, then the amplitudes a„,&_ are 
related to the amplitudes a+,b + by 



S, 



(2) 



where the 2x2 S-matrix Sr,{k,Vl) describes the trans- 
mission through the loop, thereby containing terms of 
the form e l ( kL ±<t>) ^ J%k being the longitudinal momen- 
tum of the particles, L the circumference of the loop 
and 4> the rotational phase shift. Sl is in general non- 
unitary when particles can leave the loop from another 
.junction. Time reversal symmetry (Onsager relations 
|14j ) implies that S a p = Sp a and that Sl has the form 

/ a be^ \ 
Sl = , -id, ) where a, b and c are functions of 



x be-^ c 

k. By combining eqs. |T]) and @ we find the solution 
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FIG. 1: Geometries of guided matter-wave Sagnac interfer- 
ometers: (a) MZ. Open squares mean controllable reflectivity. 
When fully reflecting we name the ifm: closed MZ. (b) 2-port 
loop. Closed squares are fully reflecting mirrors, (c) 4-port 
loop, (d) single junction loop. In all configurations horizontal 
lines mean 50-50% BSs while vertical lines mean BSs with a 
transmission amplitude it. 



where I is the 2x2 unity matrix and T = SSl, S being 
the sub-matrix of S appearing in eq. ([T]). The output 
of the ifm is obtained when propagating the amplitudes 
a + ,6 + through the ifm arms and transmitting them into 
the output port through the output junction. When no 
additional imperfections exist in the arms, the output 
amplitude is given by 



Qout 



where S' is the scattering matrix of the output junction 
and L a ,Lp are the lengths of the corresponding arms, 
such that L = L a + Lp. 

Time reversal symmetry, which determines the sym- 
metry properties of the scattering matrices S and Sl 
discussed above yields the following general form of the 
transmission through the ifm [l2j 



P(M) = Wout\ 2 = 



B + C cos(</> + (3) 



1 + D cos < 



H cos 2 



(5) 



where B,C,D,H and (3 are functions of fc. When the 
ifm has only one input and one output port, time- 
reversal symmetry and current conservation imply that 
P((j)) = P{-4>) so that f3 = or (3 = ir ("phase rigidity"). 
In what follows we show that while in an ifm having more 
than two open ports (3 may be strongly fc-dependent, such 
that an integration over a wide momentum bandwidth 
Afc washes out the cf> dependence of P(<f>), this depen- 
dence is conserved in an ifm having only 2 open ports 
and fixed f3, enabling wide momentum bandwidth oper- 
ation. 

Let us first examine a simple Sagnac ifm which does 
not satisfy the condition for " phase- rigidity" . The Mach- 
Zehnder (MZ) ifm shown in fig. la is analogous to the 



ifm implemented in [6(. It contains two 50-50% junc- 
tions where an incident particle may be either reflected 
with amplitude Sp t i — S a ^i = l/v2 or transmitted 
with amplitude S a j = Sp^i = ij\pl. This is a one- 
pass ifm, where no reflections occur from the loop ports 
a and (3 back into the loop (S — Sl = T = 0), so 
that the transmission probability at the output port o 
is P(k, 4>) = sin 2 [(0 — k5L)/2], where SL = L a — Lp is 
the length difference between the ifm arms. P(k, 4>) can 
be then written in the form ([5|) with B = — C = 1/2, 
D = H = and (3 = kSL. If SL ^ and the input flux is 
described by a Gaussian spectral distribution G(k) with 
bandwidth Afc around fc = fc, then the time-averaged 
transmission probability P(Cl) — J dkG(k)P(k, <fi(£l)) be- 
comes 



(2mA - 
1 — 77 cos — - — \l — kSL 
V h 



(6) 



where the interference visibility r\ = e~ Ak SL / 2 decreases 
with momentum bandwidth and path difference. If, in 
addition, the atomic beam is not perfectly collimated so 
that parts of it take paths surrounding different effective 
area with uncertainty SA, then rj is multiplied by a factor 
exp[-{2m/h) 2 SA 2 n 2 /2] = exp[~(SA 2 /A 2 )(/> 2 /2}. Both 
suppression factors have been observed experimentally, 
as demonstrated in fig. 3 of ref. Q. 

In contrast, we now study a closed loop ifm obtained 
by closing the ports i' and d of the MZ with mirrors 
placed in front of them (fig. la), so that each junction 
now has only 3 open ports. A particle incident from one 
of the loop arms on such a junction may be reflected 
back into the loop with amplitude = e lVij /2, where 
each phase ipij = a, (3) consist of a contribution 

9 upon reflection at the mirror and an additional 7r/2 
phase for each transmission through the BS. We then 
have S aa = -Spp = e i8 /2 and S a0 = Sp a = id /2. The 
matrix Sl has elements similar to the elements of the S- 
matrix of the output junction, {Si,)ij = S^e"^, where 
rjij describe propagation through the arms: rjij = kL ± <j> 
for i 7^ j and r/jj = 2kLj. Following the above prescrip- 
tion, we obtain the transmission probability as in ([5]) with 
(3 = 0, such that P(4>) = P(—<f>). In order to study the 
properties of the transmission when the source has a mo- 
mentum bandwidth Afc, it is insightful to realize, in view 
of eqs. ((3]) and ((U), that the output amplitude can be writ- 
ten as a sum a out = a„e^" +1 / 2 '' cL , where the am- 
plitudes a n are functions of fc with terms e lh5L , which are 
slowly oscillating relative to the fast oscillation of e mkL , 
if we assume small path length differences SL <C L. In 
the transmission probability P(k,(f>) = \a out \ 2 these fast 
oscillations will appear in cross terms a* n a n i with n =/= n' 
describing interference between trajectories with differ- 
ent number of passes through the loop in either direc- 
tion. In a realistic situation where the coherence length 
of the matter-wave source is much smaller than L (Afc 3> 
2n/L), these interference terms will be eliminated. We 
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describe this elimination by defining a slowly varying 
time-averaged transmission probability integrated over a 
period 2tt/ L of the fast oscillations 



P(k, 



k+tr/L 



k-ir/L 



dk'\a out {k', 



£M 2 > 

n=0 



(7) 

which describes the transmission of a quasi- 
monochromatic flow (Ak5L <C 1), where the coherence 
length is large enough relative to effective path length 
differences between trajectories with the same number 
of passes. In the incoherent limit where AkSL 3> 1, 
integration of P(k, 0) over the bandwidth is equivalent 
to taking the average of P over < kSL < 2ir. In 
this limit only paths with exactly the same length may 
interfere. 

For the closed MZ the transmission probability in the 
quasi-monochromatic case is found to be 



P = 



1- 1/2 (cos 2 kSL + cos 2 
1 — 1 /4(cos kSL — cos 4 



(8) 



as shown in fig. 2a for a few values of kSL (thin curves) . In 
the incoherent limit (thick curve) the visibility is not sup- 
pressed as in the simple MZ but stays fixed at rj « 0.16. 
A more general explanation for this residual visibility is 
interference between counterpropagating waves that fol- 
low trajectories with exactly equal length, as in white 
light interferometry. Change of the rotation frequency 
O changes the relative phase between these trajectories 
but they still interfere. In the closed MZ the existence of 
such pairs of interfering trajectories is allowed by inter- 
nal reflections, which in a system with a single transverse 
mode are an inevitable result of current conservation in 
a junction with only three ports. 

The rotation sensitivity of a Sagnac ifm is the minimal 
rotation frequency change 5Q m i n that generates a notice- 
able transmission change (beyond noise level). For an av- 
erage incident flux F with a Poissonian particle number 
distribution 



5Q r , 



2mA J~Ft 



dP 



(9) 



where r is the measurement integration time and S(j) min 
the dimensionless phase sensitivity per particle. The best 
sensitivity is achieved approximately near points with 
maximal derivative of the transmission. For ifm's with 
sinusoidal transmission as the MZ, the best sensitivity is 
inversely proportional to the visibility r\. 

In order to achieve better rotation sensitivity we now 
turn our attention to ifm's with a high finesse loop, where 
the finesse {(F) is defined as the ratio between the peri- 
odicity of the spectral transmission 2tt/L and the reso- 
nance bandwidth (1 — R)/\fRL, R being the probability 
to stay in the loop for a full round-trip. For a high finesse 
loop {R R* 1), J- is proportional to the average number 



of passes of a particle through the loop before exiting 
through a junction. The closed MZ ifm can be converted 
into a high finesse ifm by replacing the BS of the closed 
MZ with a "vertical" BS rotated by 90° relative to the 
MZ BS (fig. lb). A particle incident on the "vertical" BS 
has a reflection amplitude r and a transmission amplitude 
it = — r 2 , which is controllable, for example, with a 
magnetic tunneling barrier as sugg ested in hjl. T > 1 
requires that r ss 1 and t <C 1. Here the transmission am- 
plitude between the two arms S a /3 = Sp a — r, while back 
reflection through the mirror is only permitted to arm (3(3' 
with amplitude Spp = —t 2 e l6 . An analysis similar to that 
of the closed MZ yields the quasi-monochromatic trans- 
mission probability shown in fig. 2b for different values of 
kSL (thin curves). The transmission has a symmetric dip 
at (f> = 0, whose depth and width depends on the value 
of kSL. Its width for a given value of kSL is proportional 
to the spectral bandwidth A fe = (1 - r 4 )/r 2 L rj 2t 2 /rL 
of the loop, which is inversely proportional to its finesse 
T = (2ir j L) j Afc rj irr/t 2 . In the incoherent limit (thick 
curve), P{4>) rj P [l - 7]C((t))], where P rj t 2 /r, £(0) 
is a Lorenzian of FWHM ~ A k L and t] rj 0.175. 
Using eq. ©, the best sensitivity near (j) rj ±A^/2 is 
Scj) mm ~ A^/^/r, « %/2A/ry cx F^ 1 / 2 . The trans- 
mission dip at = is due to the fact that interference 
between pairs of equal length paths terminating at the 
output port is destructive when = 0. For comparison, 
we calculated the transmission of a similar high finesse 
ifm with additional ports using a combination of a 50- 
50% BS and a "vertical" BS at each junction (fig.lc). 
The quasi-monochromatic transmission of this ifm (thin 
curves in fig. 2c) is not symmetric about = and the 
visibility drops to in the incoherent limit (thick curve) . 

Finally we utilize the formalism developed here to an- 
alyze a simple ifm where each trajectory has a counter- 
propagating counterpart of the same length. The ifm 
configuration in fig. Id contains only one junction for in- 
put and output. Its output amplitude is obtained by 
substituting L a — Lp — L and S' — S in eq. Q. The 
matrix T for this ifm is diagonal, corresponding to no 
coupling between counterpropagating modes, but on the 
other hand, all trajectories of the same order n have ex- 
actly the same length, giving rise to destructive interfer- 
ence at = 0. The slowly varying transmission proba- 
bility ([7]) is then a sum over k independent contributions 
t A r 2n sin 2 [(n + 1)0] resulting in 



P(0) 



1 



COS 2 



1 + 4(r/£ 2 ) 2 sin 2 



(10) 



with a Lorenzian dip at = of FWHM A^ = t 2 /r, 
full visibility (rj = 1) and an asymptotic transmission 
P Ri t 2 /2r, as shown in fig. 2d (thick curve). The best 
sensitivity at — > is Scf> m in ~ t/{2y2r) = 
which scales as the inverse square root of the finesse 
T = 2-kt ft 2 . The transmission visibility is not affected by 
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FIG. 2: Transmissivity of Sagnac ifm's as a function of rota- 
tional phase shift <f> (a) closed MZ. (b) 2-port loop, (c) 4-port 
loop, (d) single junction loop. Thin curves in plots (a-c) 
represent the quasi monochromatic transmission P{(f)) for few 
values of phase difference between the two arms kSL, and in 
plot (d) the transmission in the presence of a scattering cen- 
ter at x = Z//3 with reflection amplitudes from r — (thick 
curve) to r = 0.25. Thick curves in (a-c) represent the trans- 
mission in the incoherent limit while all curves in (d) apply 
both in the quasi monochromatic and the incoherent limit. 
For all four models t — 0.25. As can be observed, the phase 
locking in (c) suppresses the asymmetric functions appearing 
in (b), giving rise to less smearing of the interference pattern. 
The final dip in the incoherent signal of (c) may be explained 
as interference only between exact same length trajectories, 
which due to the 7r/2 phase difference between the r and t 
amplitudes, give rise to destructive interference. 



changes in the effective path length, but may be affected 
by internal reflections from asymmetrically located im- 
perfections in the guiding potential. This affect may be 
calculated from eq. ^ when the effect of scattering cen- 
ters is included in the loop Sl matrix. As an example, we 
present in fig. 2d the incoherent transmission of the ifm 
when a single scattering center with reflectivity ampli- 
tude 0<r<0.25is placed at x = L/3. The sensitivity 
is then degraded both by visibility reduction and increase 
of the Lorenzian width, while the minimum transmission 
point stays at 4> = 0. The effect of a non-dispersive scat- 
tering center (which does not change the effective path 
length) is similar for other high finesse ifm's discussed 
above, regardless of their phase rigidity property. 

To get an estimate of the achievable sensitivity of a 
Sagnac ifm on an atom chip, we assume that an atom 
waveguide ring of radius 1cm is formed near the chip 
by magnetic field gradients of the order |VB| ~ G//_tm, 
which may be generated by wires on the chip, about 
10/jm from the surface. The centrifugal force mv 2 jr of 
the circulating atoms must not exceed /ig|VB|. This lim- 
its the maximum velocity of the atoms to v ~ lOm/sec. 
An atomic trap lifetime of about 10s permits up to 1000- 
2000 rotations, corresponding to a tunneling amplitude 



t r~j 0.035 at the BS. For the single junction loop we ob- 
tain 54>min ~ 0.013. If we assume a flux of 10 9 atoms per 
second (e.g. from a 2d MOT |l6[) we obtain a sensitivity 
of SQ ~ 5 • 10~ 43 r ad / sec/ yHz - about three orders of 
magnitude better than the best value published to date 

The Id waveguide model put forward above can be 
easily extended to a multi-mode model where a waveg- 
uide supports N > 1 transverse modes. We then have 
to replace the amplitudes a{ n , a ou t , a± and b± with N 
component vectors and the matrices S, Sl and T with 
2N x 2N matrices, which may also couple between the 
modes. The application of the suggestions in this Letter 
to real systems requires a more comprehensive study of 
the multi-mode case as well as other important effects, 
such as dispersion of ifm components, Berry's phase in 
a magnetic ring potential [l8| and atom-atom collisions. 
The ideas raised in this work concerning time reversal 
symmetry in a rotating quantum system with many pos- 
sible paths are closely related to the theory of weak local- 
ization and coherent back scattering in mesoscopic elec- 
tronic systems in the presence of a magnetic field [Toj]. 
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